ON MOMENTS OF THE PARABOLIC ANDERSON MODEL 
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Abstract. We study the parabolic Anderson model with nearest neighbor jumps and space-time 
white noise (discrete space / continuous time). We prove a contour integral formula for the second 
moment and compute the second moment Lyapunov exponent. For the model with only jumps to 
the right, we prove a contour integral formula for all moments and compute moment Lyapunov 
exponents of all orders. 



1. Introduction and main results 

1.1. Nearest-neighbor parabolic Anderson model. The nearest-neighbor parabolic Anderson 
model on Z is the solution to a coupled system of diffusions on [0, oo) given by 

j t Z p {t, n) = A^Zpit, n) + pZp(t, n)dW n (t). (1) 

Here t G M+, n € Z, and the operator A p ' g (which is the generator for a nearest neighbor continuous 
time random walk) acts on functions of n as 

A™f(n) = Pf(n - 1) + qf(n + 1) - (p + q)f(n). (2) 

We will assume that p, q > and p + q = 2. The collection {W n (-)} n< =z are independent Brownian 
motions and (3 £ R+. For a general physical and mathematical background on this model, the 
reader is referred to [5j [9] as well as part I of [8] for more recent developments. 

The coupled diffusions can be considered as modeling population growth in a random, quickly 
changing environment at each spatial location, and with migration between locations. Consider a 
population of many small particles living on the sites of Z. There are three forces acting upon this 
system: 

(1) Each particle at time t and lattice site n independently duplicates itself at rate r + (t,n); 

(2) Each particle at time t and lattice site n independently dies at rate r_(t,n); 

(3) Each particle at time t and lattice site n independently jumps to a neighboring site n — 1 
with rate q and n + 1 with rate p. 

Letting m(t, n) be the expected population size at time t and location n, one finds that [5] 

— m(t, n) = A p ' q m(t, n) + (r + (t, n) — r_(t, n))m(t, n). 

If the duplication and death rates are independent in space and quickly mixing in time, the factor 
(r+(i, n) — r_(t, n)) is well modeled by j3dW n (t) where j3 modulates the relative rates of jumping 
and duplication/death. 

Closely related to the above branching diffusion representation, the Feynman-Kac representation 
for this coupled system of diffusions writes Za(t, n) as point to point partition functions for a random 
polymer model (for background on these models, see the review [6]): 



Zp(t,n) = £. 



7r(0)=0 



(3) 



i 
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where 7r(s) is a Markov process with state space Z and generator given by A q ' p (which is the adjoint 
of A p ' q ), and £' 7r (o)=o is the expectation with respect to starting 7r(0) = 0. We write E for the 
expectation over the disorder. 

1.2. Lyapunov exponents and intermittency. We consider two variants of the Lyapunov expo- 
nents for the parabolic Anderson model. Consider a velocity v € K. Then the almost sure Lyapunov 
exponent with respect to velocity v is given by 

71 iP\ = lim 7 lo g^/s(*, NJ)- 

t— >oo t 

The existence of this almost sure limit is due to a sub-additivity argument (see [5], Section IV. 1). 
The p-th moment Lyapunov exponent with respect to velocity v is given by 



lim — log E 

t— >oo t 



(Zp(t, [ut\)Y 



If the initial data Zp(Q,ri) is stationary with respect to shifts in n, then the exponents are, in 
fact, independent of v. We, however, will consider initial data in which Zp(0,n) = l n =o and hence 
the exponents will depend on the velocity v non-trivially. 

Intermittency is defined as the strict ordering of exponents: 



7i (p; i/) < 7i (p; < — ~ — < 



< 



2 3 

The weak ordering of exponents is a consequence of Jensen's inequality (for the first inequality) and 
Holder's inequality (for all subsequent inequalities). A useful fact is recorded in the following (cf. 
[5], Theorem III.1.2): 



Lemma 1.1. If for any k > 1, 



then for all p > k 



7fc(P» < 7fe+i(/3;^) 



k 



< 



P 



k + 1 
p+1 



(4) 
(5) 



Proof. Holder's inequality implies that 

^ 2 



E 



(z p (t, [ut])y 



< E 



(Zp(t, [vt\)) 



k+h 



E 



f Mt, [vt\)) 



k-h 



from which follows the inequality 



If @ holds, then 



k 



(6) 



"1k{P\v) < ~i ~ _ r7fc+i(/3; v) 



k + 1 

Plug this into the right-hand side ([6]) with h = 1 and k replaced by k + 1. Rearranging terms yields 
([5]) for p = k + 1. This can be repeated inductively, and yields the claimed result. □ 



Intermittency is related to the presence of high peaks of Zp(t,n), with low probability. Fix a 
zh that 7fc( f ; ^ < a < 7fc +^ ;I/) ; then we know that W(Zp(t, vt) > e at ) > as t -+ oo. Writing 



E 



Z (t, [ut\) 



k+1 



E 



Z /3(t, Y u \ t ) k+ll Z fi (t,\ut\)<e a 



+ E 



Z p (t, [vt\) k+l l Zp[t ^ )>t 
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we observe that the first term is < e a ( k + 1 ) t ^ but the sum of the two is asymptotically e yk+1 ^' v ' t , 
which is exponentially (as t grows) larger than e a ( k + 1 ) t . This means that the event {Zp{t, vt) > e at } 
gives overwhelming contribution to the k + 1 moment. On the other hand, 



E 

hence 



Z (t, [vt\) k ] > e akt ¥ (Z p (t, [ut\) > e at ) 



(Zp(t, L^j)>e^)<^^=exp(-fa 



gafct IV k 

which is exponentially small. 

In the case of spatially translation invariant ergodic solutions Zg(£,n), the consequences of in- 
termittency may be interpreted via spatial averages over large balls at a fixed (large) time. Thus, 

one can talk about islands where the solution is at least e (as opposed to the typical value of 

e i\(p,v)t^ w bose spatial density is not more than e \ fc+1 fc / . Our results are for delta initial 
data Zj3(0,n) = l n =o and not stationary initial data. However, Lyapunov exponents are remarkably 
robust with respect to different initial data ([5], Lemma III. 1.1) and one can hope, however, that 
the information on 7/c(/3;z/), for stationary solutions can be read off from non-stationary ones as 
well. 

1.3. Main results. All of our results pertain to the nearest-neighbor parabolic Anderson model 
with delta initial data: Z^(0,ra) = l n= o- Our first result is a formula for the two-point moment of 
the model. 



Theorem 1.2. For n\ > n 2 , 

2 

~[Zp(t,rii 



E 
where 



i=l 



1 11 {PZ1 ~ ffV {PZ2 -jf\ M F%M)WM--, (7) 



(2tt0 2 J J (pz t - qz^ 1 ) - (pz 2 - qz 2 x ) - 2/3 2 hni " ' £ '™ 2V ' z x z 2 



and where the contour of z\ is the unit circle and the contour for z 2 is a circle around of radius 
sufficiently small so as not to include any poles of the integrand aside from z 2 = 0. 

This theorem is proved in Section [2j We state an immediate corollary of this theorem which 
shows how it yields an exact solution to a renewal equation. Let p(t, n) denote the heat kernel for 
the generator A p ' q (i.e., the solution to ^p(t,n) = A p ' q p(t,n) with p(0, n) = l n= o). 

Corollary 1.3. For all t G R + , n£Z, and (3 G M + the equation 



f(t,n) 2 =p(t,n) 2 +(3 2 / ds p(t — s,n — m) 2 f(s,m) 2 

J° m=-oo 

with initial data /(0,n) = l n =o is solved by the right-hand side of ([?J) with n\ = n 2 = n. 

Proof. The integral formulation of (fT]) is 

ft °° 

Zp(t,n) =p(t,ri)+ ^ p(t - s,n - m)Z l3 (s,m)dW m (s). 

■'° m=-oo 

Square both sides and take expectations. Call f(t, n) = E [Zp(t, n) 2 ] , then the ltd isometry implies 
the claimed result. □ 
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Via asymptotic analysis, Theorem 11.21 enables us to calculate the first and second moment Lya- 
punov exponents. 

Theorem 1.4. The first moment Lyapunov exponent at velocity p — q, for the nearest-neighbor 
parabolic Anderson model is given by 7i(/3;p — q) = 0. The second moment Lyapunov exponent at 
velocity p — q, for the nearest-neighbor parabolic Anderson model is given by 

72(/3;p - q) = H(z ) 

where 

H{z) = -z{ps{z) + q(s(z)y 1 - 2 - (p - q) log (s(z)) + pz + qz~ l - 2 - (p - g) log z 

with 



(pz- qz~ l + 2/3 2 ) + J(pz - qz" 1 + 2/3 2 ) 2 + 4pq 

S{Z) = 2~p 

and where zq is the unique solution to H'(z) = over z G (0, oo). 
When p = q = 1, 

z = \ (-/3 2 + V4 + /3 4 ) , s(z ) = \ [f + ^4 + /3 4 ) , H(z ) = 2(^4 + ^- 2), 
which implies that 

72 (/5;0) = 2(V4 + /3 4 -2) 
/or t/ie standard (p = q) parabolic Anderson model. 

This theorem is proved in Section [2] via asymptotic analysis of Theorem 11.21 We include the full 
details only for the case p = q = 1. 

Remark 1.5. The above theorem is stated only for a velocity given by v = p — q. This is because 
in that case, the first Lyapunov exponent is 0. For general p — q ^ the same approach as given in 
the proof provides the exact values of the first and second moment Lyapunov exponents. 

We now turn attention to the one-sided case of the nearest-neighbor parabolic Anderson model, 
where p = 2 and q = 0. In this case we may extend the result of Theorem 11.21 to arbitrary joint 
moments. For k > 1, define 

W> = {n= (n 1 ,n 2 ,...,n k ) G (Z >0 ) k : m > n 2 > . . . > n k > 0}. (8) 



Theorem 1.6. For all k > 1 and ft G W> , 



E 



~[Zp(t,ni) 



i=l 



n z a - z b -A- e*( Zl l ) dzi 
* _ y, _ R2 11 



(27u) k J J 11 z a - z b - ft 2 J- 1 4 

V ' J J l<a<b<k a H i=l * 



where the integration contour for z a is a closed curve containing and the image under addition by 
f3 2 of the integration contours for 2% for all b > a (for an illustration of possible contours see Figure 

EP. * 

This theorem is proved in Section [3l Asymptotics of this formula yields all the moment Lyanpunov 
exponents. By Brownian scaling it suffices to consider just /3 = 1. 

Theorem 1.7. For any k > 1 and v > 0, the k-th moment Lyapunov exponent at velocity v for the 
one-sided (p = 2 and q = 0) nearest-neighbor parabolic Anderson model with /3 = 1 is given by 

7fc(l,i/) = H(zq) 
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Figure 1. Valid contours for equation (|14p with k 
the z\ contour contains the image of Z2 plus /3 2 . 



2. The inner contour is Z2 and 



where 



fk-\ 



\i=0 / 

and where zq is the unique solution to H'{zq) = with zq £ (0,oo). 

This theorem is proved in Section [3] via asymptotics of Theorem 1 1.61 Figure [2] records the plot of 
the various Lyapunov exponents. 

The almost sure Lyanpunov exponent was conjectured in [13] and proved in [12] . and it is given 
by (for /3 = 1) 

7i(l;i/) = -- + inf(t-i/tf(t)) 

where *&(t) := [logr]'(t) is the digamma function. 

There are two ideas which are behind the results of this paper. The first idea is the content of 
Propositions 12.2 1 and 13. ll which show that one can compute the moments of the parabolic Anderson 
model via solving a system of coupled ODEs with spatial variables n £ W> and specific boundary 
conditions. This reduction to solving ODEs on W> only works for k = 1,2 with the general p, q 
nearest-neighbor model. However, for p = 2 and q = 0, the reduction holds for all k. The second 
idea is that the system of ODEs can be explicitly solved via a certain nested-contour integral ansatz 
that originated from [3]. This is the content of Propositions 12.3 1 and 13.31 

The rest of the paper is as follows: In Section [2] we show how the moments of the parabolic 
Anderson model can be computed via a coupled system of ODEs. We then solve this system 
and use this solution to prove Theorems 11.21 and 11.41 In Section [3] we show how in the one-sided 
model, all moments can be computed via ODEs and we provide integral formulas which solve these 
ODEs. From this we are able to prove Theorems 11.61 and 11.71 We conclude the section with a 
non-rigorous replica trick calculation (used extensively in the physics literature) and show how from 
this calculation one recovers the almost sure Lyapunov exponent for this model. In Section [4] we 
briefly discuss the continuous space parabolic Anderson model (i.e., the stochastic heat equation 
with multiplicative noise) and record its moments and Lyapunov exponents. 
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1.5- 




-0.5 J 



FIGURE 2. Plot of one-sided parabolic Anderson model Lyapunov exponents ver- 
sus the velocity v (which plays the role of the diffusion constant). Here we have 
normalized Zi(t, [_vt\) by the zero noise solution, so that 7i(l;z/) = 0. The lowest 
curve in the plot is the normalized 71 (1; u) and the higher curves are the normalized 
7fe(l;i/)/A; (increasing in height with k). This demonstrates the intermittency of this 
parabolic Anderson model and the shape of this plot is very similar to those on page 
105 of 0. 

1.4. Acknowledgements. We thank Frank den Hollander for discussions on the parabolic An- 
derson model and useful comments on an early draft of this article; and Davar Koshnevisan for 
pointing out Corollary 11.31 AB was partially supported by the NSF grant DMS-1056390. IC was 
partially supported by the NSF through PIRE grant OISE-07-30136 and DMS-1208998 as well as 
by Microsoft Research through the Schramm Memorial Fellowship, and by the Clay Mathematics 
Institute through the Clay Research Fellowship. 

2. Nearest neighbor parabolic Anderson model 

The first step in our computation of the moment Lyapunov exponents of the parabolic Anderson 
model is the following reduction to a coupled system of ODEs with two-body delta interaction. 
Recall the definition of the nearest-neighbor parabolic Anderson model Zp(t,n) and the operator 
A p ' q given in the introduction. Write [A p,9 ] i for the operator which acts as A p ' q on the i-th spatial 
coordinate. 

Proposition 2.1. Asumme v : R+ x Z fc — > R solves 
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(1) For all n G Z fc and t G R + , 

j t v(t;n) = Hv(t;n), H = § £ [A**], + ±/? 2 £ l na=n6 ; 

j=l a, 6=1 

(2) For aZZ permutations of indices a G Sk, v(t; an) = v(t; n); 

(3) For all n G Z fc , lim t ^ n) = U i=1 1^=0/ 

(4) For all T > 0, i/iere exists c, C > sitc/i i/iai /or a// n G Z fc em<i a// i G [0, T], 

\v{t;n)\ < ce c W n W\ 
Then for n G Z k , v(t; n) = E J]Li "0 • 

Proof. This result is well-known and can be found, for instance, in Proposition 6.1.3 of [3]. The 
purpose of the fourth hypothesis on v is to ensure uniqueness of solutions to the system of ODEs 
given by the first three hypotheses. The fact that this exponential growth hypothesis is sufficient 
for uniqueness can be proved in the same manner as given in the proof of Proposition 4.9 in 

One way to see why this should be true is to consider the Feynman-Kac representation for Zg(t, n) 
which is given in equation ([3]). The k factors of Zp lead to k paths. The expectation E over the 
Gaussian disorder (white-noise) can be taken inside the path expectations £ and calculated exactly 
yielding the exponential of the pair-wise local time for the k paths. This accounts for the delta 
interaction seen above. □ 

It is a priori not clear how one would start to solve the system of ODEs in the above proposition, 
one reason being that it is inhomogeneous in space. An idea from integrable systems (related to 
the coordinate Bethe Ansatz) is to instead try to solve a homogeneous system of ODEs and put the 
inhomogeneity into a boundary condition. If the number of boundary conditions is k — 1, then there 
is generally hope in solving the system by combining fundamental solutions of the homogeneous 
system in such a way that the initial data and boundary conditions are met. 

For the general p, q case, it appears that this reduction to k — 1 boundary conditions only works 
when k = 2 (in which case there is just one boundary condition). When p = 2 and q = the 
reduction works for all k (see Section [3]). 

Proposition 2.2. Assume u : R + xZ 2 ->l solves 



(1) For all n G I? and t G 



i=l 

(2) For n such that n\ = 712 = n 

P Q P Q 

Tpu(t; n) := (3 u(t; n, n) + -^ u {t', n,n — 1) + - u(t; n + 1, n) — -^u{t\ n — 1, n) — -^ u (t\ n > n + 1) = 0; 

(3) For all n G Z k such that n\ > 122, lim^o u(t; n) = Yt£=l l«i=0;' 

(4) For all T > 0, there exists c, C > such that for all n G Z fc such that n\ > n-i and all 
t€ [0,T], 

\u(t;n)\ < ce c ^\ 
Then for neZ 2 such that n\ > ri2, u(t; ft) = v(t; n) = E H?=i ^p{^i n i 
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Proof. We show that restricted to n \ >ri2, u(t;n) symmetrically extended to Z 2 solves the system 
of ODEs in Proposition 12.11 and hence u(t;n) = v(t;n). For n\ > ni it is clear that u and v solve 
the same equation. For n\ = ni = n, 



dt' 

= (/3 2 — 2)u(t; n, n) + pu(t; n, n — 1) + qu(t; n + 1, n), 



-u(t; n,n) = | (pu(t; n — 1, n) + quit] n + 1, n) + pu(t; n, n — 1) + qu(t; n, n + 1) — 4u(i; n, re) J 



where the second line followed from the relation imposed by the assumption (2). Now compare this 
to the equation v(t; n, n) satisfies: 

—v(t;n,n) = i (pv(t; n — 1, n) + qv(t; n + 1, n) — 2v(t; n, n)^ 

+ \ (p v (t'i n,n — I) + qv(t; n, n + 1) — 2v(f; re, n)J + (3 2 v(t; n, n) 

= (/3 2 — 2)i>(t; n, n) + yw(i; n, n — 1) + gi>(f; n + 1, n), 

where the second line followed from the symmetry hypothesis on v. Observe that on the diagonal 
n i = n>2, both u and v solve the same equation. Therefore they solve the same equation for all 
n\ > n2 and hence (since the other hypotheses of Proposition 12.11 are satisfied) u(t; n) = v(t; ft). □ 

We may now explicitly solve the system of ODEs defined in Proposition 12.21 

Proposition 2.3. For k = 2 and n\ > 112, the system of ODEs in Proposition \2.2\ is uniquely solved 
by 

u(t;n lina ) = J^ff {pzi _ _ {pz2 _ _ 2 p *t,nM)F t ,r»(*2)--> (9) 



where 



F™{z) = z -« e |(P^- 1 -2) 



t,n 

and where the contour of Z\ is the unit circle and the contour for Z2 is a circle around of radius 
sufficiently small so as not to include any poles of the integrand aside from Z2 = 0. 

Remark 2.4. The right-hand side of ([9]) is easy to generalize to all k (see e.g. Section 6.1.2 of [3] 
or Proposition 13.31 below in the totally asymmetric case where p = 2 and q = 0). However, it is not 

at all clear if such a generalization would have anything to do with E Y\i=i Zp{ti n i 

Before proving this proposition, we note that Theorem 11.21 follows as an immediate corollary of 
the above result and Proposition! 



Proof of Proposition \2.3l We prove this proposition by checking the hypotheses of Proposition 
Hypothesis 1 follows from the Leibniz rule and the observation that 

= A»«i$?(*). 

To check hypothesis 2, we apply Tp to u(t; n) when n\ = n2 = n. The operator Tp can be taken 
inside the integration. It acts on Ff^{z{)Ff^{z2) as 

Tp(F™( Zl )F™(z 2 )) = -lF^(z 1 )F^(z 2 )((pz 1 - qz^) - (pz 2 - qz^) - 20< 
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The factor ( (pz\ — qz 1 1 ) — (pz 2 — qz 2 1 ) — 2/3 2 j cancels with the same term in the denominator of 
the integrand, yielding 

Tpu(t,n) = -J^pfffa - qz^) - (pz 2 - qz^))F™( Zl )F™(z 2 )^^-. 

Since the integration contours are the same for both Z\ and z 2 and since {{j>z\ — qz^ 1 ) — (pz 2 — qz^ 1 )^ 

is antisymmetric in z\ and z 2 , while the rest of the integrand is symmetric, one immediately sees 
that the integral is zero as desired to check hypothesis 2. 

Hypothesis 3 is checked via residue calculus. The t — > limit can be taken inside the integrand 
and we are left to show that for m > n 2 , 

1 'c( Zl ,z 2 )z^z^^^ = l„ 1=0 l„ 2 =o, (10) 



(2^) 2 J J 1 1 z 1 z 2 

where 

r , \ (p«i - i z i l ) - (p z 2 - qz^ 1 



(pzi - qz 1 *) - (pz 2 - qz 2 *) - 2/3 2 ' 
and where the contour of z\ is the unit circle and the contour for z 2 is a circle around of radius 
sufficiently small so as not to include any poles of the integrand aside from z 2 = 0. 

(1) If n 2 < then in (I10p we may shrink the z 2 contour to 0. Observe that for z\ fixed on the 
specified contour, C(z±, z 2 ) is analytic in z 2 in a small neighborhood of z 2 = 0, with a value 
of C(zi, 0) = 1. The term z^" n2 ^- does not have a pole at (because n 2 < 0) and hence in 
this case u(0; n) = 0. 

(2) If n 2 = then in (|10p let us shrink the z 2 contour to 0. The term z^™ 2 ^ has a simple pole 
at and hence the integral evaluates as 

1 I -n! dz l , 

(3) If n 2 > this implies that n\ > as well. Then we can expand z\ to infinity. As we do this, 
we encounter a pole at z\ such that 

{pzi - qz^ 1 ) - (pz 2 - qz^ 1 ) - 2(3 2 = 0. 

For each z 2 there is only one such pole r(z 2 ) which comes when z\ ~ p 2 ^ 1 l° r sma ll \ z 2\- 
The reason why only one pole is crossed is because the other pole coming from this term 
is approximately z 2 , which is already contained inside the Z\ contour. Before analyzing the 
residue, observe that because m > 0, there is at least quadratic decay in z\ at infinity, so 
there is no pole at infinity. Thus, the integral in z\ is given by its negative residue at r(z 2 ). 
The negative residue at Z\ = r(z 2 ) is evaluated as 

1 f 2f3 2 r(z 2 ) ni _ x i 



v iz 2 )Y n ^z- n ^dz 2 (11) 
2-kl J pr(z 2 ) + q(r{z 2 )) iV 

where the integral in z 2 is over a small circle around the origin. It is easy to see that 

W 2 r(z 2 ) 

pr(z 2 ) + g(r(z 2 )) -1 

is analytic in a neighborhood of z 2 = and its value at z 2 = is 2/3 2 /p. Thus the entire 
integral (fTTI) can be evaluated as the residue at z 2 = 0. Since n\ > n 2 , there is, in fact, no 
pole at z 2 = 0, thus the integral equals 0. 
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Combining the above cases we see that the only case in which u(0; ft) is non-zero is when n\ = ri2 = 0, 
in which case it is 1. This confirms the initial data of hypothesis 3. 

Hypothesis 4 follows via easy bounds of the integrand of u(t; n). □ 

Having proved the two-point moment formula for the nearest-neighbor parabolic Anderson model, 
we can now extract the second moment Lyapunov exponent via asymptotic analysis. 



Proof of Theorem \l-4\ We will present a complete proof only in the case p = q = 1 since this 
simplifies the (rather technical) analysis. For the moment we keep the p and q and only set them 
equal when necessary. 

Let us start by proving 71 = from the formula 



which one easily checks via either Proposition 12.11 or I2.2L Let n = |_(p — q)t\ and observe that (up 
to an insignificant correction coming from the fractional difference between n and (p — q)t) 

1 f dz 

E[Zp(t,n)] = — <p e tG{z) — , G(z) = pz + qz~ l - 2 - (p - q) log z. 

ZTTt J Z 

We want to study this as t — > 00, thus we can use the standard Laplace's method (see Lemma 1 5. II 
for t = 1) to perform the asymptotics. The critical point equation for G(z) is 

G'(z) = p — qz~ 2 — (p — q)z~ 1 = 

which is solved by z = 1 or z = —q/p. The critical point z = 1 corresponds to the larger value 
of G(z), namely G(l) = 0. Observe that we can deform the z contour to lie on the unit circle e id . 
As a function of 9 along this contour Re[G(z(0))] = 2cos(#) — 2. This shows that the real part of 
G(z) decays monotonically with respect to the angle 6 away from z = 1. In the vicinity of z = 1, 
Re[G(z)] decays quadratically in the imaginary directions. Invoking Lemma 15.11 for £ = 1 shows 
that 

7l := lim - log (E [Z fi (t, (p - q)t))) = 0. 

t— >oo t 

To calculate 72 we use the formula for E [Zp(t,ni)Zp(t,ri2)] proved in Theorem 11.21 In order to 
perform the asymptotics we would like to deform the Z2 contour to coincide with the Z\ contour. 
While doing this we encounter a pole and must take into account the associated residue, in addition 
to the evaluation of the remaining integral on the new contours. The result of this manipulation is 

E [Zp(t, (p - q)tf] = (A) + (B) 

( a\ - 1 II (pzi ~ qzY 1 ) - (pz 2 - qz^ 1 ) m p>q dzi dz 2 

with the Z2 contour coinciding with the z\ contour, and 

(B) = — t fr^i- m{Z2) —- (12) 



where 



2ttl J p + q[s(z2)] 2 z 2 

Note that this residue term comes from Z\ = 5(^2) where s(z 2 ) is given in the statement of the 
theorem. Since the definition of s(z) involves a square-root, for z complex we specify that for 
z = re Ld with 9 £ (-00,00), \fz = yjre id l 2 . 
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As it will turn out, the residue term (B) has a larger exponential growth rate than the integral 
term (A) and hence accounts entirely for the value of 72. To see this, we compute the exponent for 
both terms. We claim that 

lim - log [(A)] = 0. 
t— yoo t 

It is easy to see why this is true. The contours for z\ and z 2 can be chosen so that there exists a 
positive constant C such that along the contours 

\9(zi,z 2 )\ < C 

where 

(pzi - qz^ 1 ) - {pz 2 - qz 2 x ) 

g(zi,z 2 ) = — — — — 7^- 

[pzi -qz 1 ) - (pz 2 -qz 2 ) - 2f3 2 
Likewise, in a neighborhood of (zi,z 2 ) = (1,1), one checks that 

c\z\z 2 - z 2 - z\z\ + zi\ < \g(zi,z 2 )\ 

for a small, yet positive c. One easily checks the remaining assumptions necessary to apply Lemma 
15.11 (with £ = 2) and therefore finds that as t — > 00, the growth of the integral defining (A) is 
governed by the value of G{z\) + G(z 2 ) at the critical point (1, 1). By comparison to the calculation 
performed above for 71, we find that 

lim ilogp)]=2 7 i. 

t— >oo t 

Since 71 = 0, the claimed result for (A) follows. 

Turn now to the residue term (B) and call z 2 simply z. From this point on we will assume that 
p = q = 1 to simplify the analysis. A similar, albeit lengthy, analysis can be performed for all p and 
q. Notice that 

^0 = l 2 (-/? 2 + V^+P 1 ) 

is a critical point of H (z) and that the contour for z can be deformed without crossing any singu- 
larities of (fT2|) to the contour F parameterized by z = z§e id for 9 G [0, 2tt]. We wish to use Laplace's 
method by applying Lemma 15. II with £ = 1. It is straightforward to check hypotheses (1), (3) and 
(4) of the lemma. Hypothesis (2) requires more work. 

To check hypothesis (2a) of Lemma [5J] observe that H'(zq) = 0, H"(zq) 7^ and H(z) is analytic 
in a neighborhood of zq. Thus it immediately follows that it behaves locally like H(zq) + c(z — 
zq) 2 + o((z — zq) 2 ). Hypothesis (2b) requires that 

p{6) := 2Re[H(z) - H(z )} 

(the factor of 2 is irrelevant) is strictly negative for all z S T \ {zq}. In fact, by symmetry of H 
through the real axis, p{9) = p(2n — 9) and hence this strict negativity needs only be checked for 
9 G (0, 7r]. By utilizing the fact that 



we find that 



p{9) = V / 4 + /3 4 (cos(0)-2)+Re 
Since for a, b real, 



\l /3 4 (2 - cos(0)) 2 - (4 + /3 4 ) sin(#) 2 + 4 + i2/3 2 (2 - cos(0))\/4 + /3 4 "sin(#) 
ReV a + ib = \ - — — ^ , 
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checking hypothesis (2b) reduces to checking that for 6 G (0, n] 



p(6) = v/4T^(cos(0) - 2) + y < 

where 

a = /3 4 (2 - cos(#)) 2 - (4 + /3 4 ) sin(0) 2 + 4, b = 2/3 2 (2 - cos(0)) + £ 4 sin(0). 
This inequality is equivalent (via elementary algebraic manipulations) to showing that 

-16(4 + /3 4 ) ( - 36 - /3 4 + 32 cos(#) + (-4 + /3 4 ) cos(20)) sin(#/2) 2 > 0. 

Letting x = cos(6) the above inequality is equivalent to showing that 

f(x) = -2(l-x 2 )/3 4 (4 + /3 4 )- (4+(-2 + x) 2 /3 4 -(l-x 2 )(4 + /3 4 )) 2 (13) 

+ (4 + (-2 + x) 2 (3 A - 2(2 - x) 2 (4 + /3 4 ) - (1 - x 2 )(4 + /3 4 )) * > 

for x G [—1, 1). The above equation is of fifth degree in x, however one checks that x = 1 is a 
root. Direct inspection shows that all four other roots are off the real axis (in fact they can be 
computed explicitly since after factoring (1 — x) the equation is quintic). Thus, to show the desired 
inequality, it suffices to show that the derivative of the right-hand side of (|13p at x = 1 is negative. 
Computing, one sees that /'(l) = —4(64 + 20/3 4 + /3 s ) which is clearly negative. Having shown this 
final inequality, it immediately follows that p{6) < for all G (0, n], just as necessary to prove 
hypothesis (2b) of Lemma 15. 11 

Applying Lemma ED to (|12h we find that 

lim jlog[(A)]=H(z ). 

Since H(zq) is positive (as compared to the contribution from the (A) term asymptotics) we conclude 
that 7 2 (/3;0) = H(z ). □ 



+ Va 2 + b 2 



3. The one-sided parabolic Anderson model 

■ 



We now focus on the one-sided case where p = 2 and q = 0. Recall the definition of W> given 



in (TST) and the definition of A p ' q given in ©. In particular, A 2 >°/(n) = f(n - 1) - f(n) 

Proposition 3.1. Fix k G Z>o- 

Part 1. Assume v : R+ x (Z>_i) fe solves: 



(1) For all n G (Z>_i) fc and t G 



^v(t;n) = Uv(r;n), H = \ £ [A 2 <°] 4 + I/3 2 £ l na=n ,; 

j=l a, 6=1 



(2) For all 'permutations of indices a G Sk, v(t; an) = v(t; ft); 

(3) For all n G W$ , v(0;n) = Y\ti 1^=0/ 

(4) For aZZ ?i G (Z>_i) fc such that rifc = -1, n) = for all t G E+. 
T/ien /or a// n G W^ 0J E rji=i n = ™)- 

Part Assume u : M. + x (Z>_i) fc — >■ M solves: 
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(1) For all n £ (Z>-i) k and t £ R + , 



k 

ju{t-n) = \^[^].u{t-n)- : 

i=l 



(2) For all n £ (Z>_i) fc such that for some i £ {1, . . . , k — 1}, ?ij = n^i, 

([A 2 '°]. t -[A 2 ' ] m -2^)n(t;n) = 0; 

(3) For all ft £ (Z>_i) fc suc/i £/iai = — 1, u(t;ft) = /or a// £ £ K + ; 

(4) For all n £ W* , u(0; n) = ft ti ^=0- 



Then for all n £ W> , E ElLl Z /9(*> n ») = u(t;n). 
Proof. This is contained in Proposition 6.3 of [1]. □ 

Remark 3.2. Part 1 of the proposition is essentially a specialization of Proposition 12. II to the case 
p = 2, q = 0. However, due to the delta initial data and the one-sided nature of the operator 
A 2 ' , v(t;n) with ft £ { — 1, 0, . . . , m} k evolves autonomously as a closed system of ODEs. This 
ensures uniqueness of solutions and explains why we no longer require the at-most exponential 
growth hypothesis which was present in Proposition 12.11 Part 2 of the proposition is an extension 
of Proposition 12.21 to general k, but only for p = 2, q = 0. The fact that this holds for all k is what 
enables us to solve for higher than second moments in this one-sided case. 

The system of ODEs in part 2 of Proposition 13.11 can be solved via a nested-contour integral 
ansatz introduced in [3] and further developed in [3]. This yields the following generalization of 
Proposition 12.31 to all k but only for the one-side (p = 2, q = 0) case. 

Proposition 3.3. For all k > 1 and n £ W> , the system of ODEs in part 2 of Proposition \3.1\ is 
uniquely solved by 

^ > J J \< a <b<k f i=l * 

where the integration contour for z a is a closed curve containing and the image under addition by 
p 2 of the integration contours for Zb for all b > a. 

Before proving this proposition, we note that Theorem 11.61 follows as an immediate corollary of 
the above result and Proposition 13.11 Part 2. 

Proof of Proposition \3.3[ We check the four conditions for the system of ODEs in part 2 of Propo- 
sition [3J3 Condition (1) follows by Leibnitz rule and the fact that 

d e t{ - z ~^ _ 1 e'^" 1 ) 
dt z n ~ 2 z n 

for all z £ C \ {0}. 

Condition (2) follows by applying ^[A 2 ' ]^ — [A 2 '°] i+1 — 2/3 2 ^ to the integrand of the right-hand 

side of ()14p . The effect of this operator is to bring out an extra factor of 1{z% — Zi+\ — [3 2 ). This 
factor cancels the corresponding term in the denominator of the product over a < b. Without the 
pole associated with this term, it is possible to deform the Zi and Zi+\ contours to coincide and 
since rij = rij+i, we find that 

([A 2 '°] t - [A 2 - ] . +1 - 2/3 2 ) u(t;n) = I d Zi I dz t+l {z t - z i+1 ) f \z t ) f (z i+1 ) , 
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where f(z) includes all of the other integrations aside from those in Zi and The above integral 
is clearly by skew-symmetry, thus confirming condition (2) as desired. 

Condition (3) follows by observing that if = —1 then in the right-hand side of (|14p there is 
no pole at in the Zf. variable. By Cauchy's theorem, this means that since the Zf~ contour only 
contains and no other poles of the integrand, the entire integral is 0, as desired. 

Condition (4) follows from three easy residue calculations. From above, if rik < 0, the integral 
in (I14p is zero. Similarly, if n\ > the integrand in (|14p has no pole at infinity and the z\ contour 
can be freely deformed to infinity. By Cauchy's theorem this means that the entire integral is zero. 
The only possible nonzero value of u(0; n) is (due to the ordering of the elements in n £ W> ) when 
n\ = ■ ■ ■ = Uk = 0. The value of u for this choice of n is readily calculated via residues to equal 1, 
just as desired. □ 

We now show how asymptotics of the result of Theorem 1 1 . 61 yields a proof of the moment Lyapunov 
exponents claimed in Theorem 11.71 



Proof of Theorem \l. r /\ The starting point of this proof is the moment formula of Theorem 11.6 
Setting all rii = \yt\ and /3 = 1 we find that 



E 



l[Zp(t, \yt\) 



(2ttlY 



n j^rii < 15 > 

l<a<b<k i=l 



where the integration contour for z a is a closed curve containing and the image under addition by 
1 of the integration contours for Zf, for all b > a. From now on we will study the right-hand side of 
the above equality, with F^'P t , (z^ replaced by F t 2 ^(zj), as the asymptotic effect of this modification 
is easily seen to be inconsequential. 

In order to perform the asymptotic analysis necessary to compute the moment Lyapunov expo- 
nents, we would like to deform our contours to all coincide so as to apply Lemma [5. 11 This requires 

2 

deforming all contours to pass through a specific critical point of logF t ' ut (z). However, due to the 
nesting of the contours, such a deformation requires passing a number of poles. The following lemma 
records the effect of such a deformation. 

Lemma 3.4. Consider a function f(z) which is meromorphic and has only one pole, located at 0. 
For k > 1, set 

1 I In za- zb /Oi) • • • f{zk) 



fJ-k 



dzi ■■■dz k 

, , . ZA — Zu — 1 Z\ ■ ■ ■ Zh 

J \<A<B<k A D 



where the integration contour for z a contains and the image under addition by 1 of the integration 
contours for Zb for all b > a. Then 



A=1 m 12 m 2 ... hfc 



where 



/A = mi !m 2 !...(2^) ^'"^ det 



Xi + Wi — Wj 



-|f(A) 

[] f( Wj )f( Wj + !)■■■ f( Wj + Xj - l)dwj. 

i,j=l j=l 

(16) 

Here X = l m i2 m2 • • • h k means X is a partition X = (Ai > A 2 > • • • > 0) such that ^ A, = k; mi 
records the number of entries of X equal to i; £(X) is the number of non-zero entries in X; and for 
1 < j < £{X) the Wj contours are all chosen to be the same contour as zj^. 
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Proof. This is given in [3] as Proposition 6.2.7. □ 

The following lemma will also be helpful in completing our asymptotics. 

Lemma 3.5. Consider a (finite) collection of functions I\(t) with A G A and |A| < oo. Assume 
that there exist 

7 A = lim -log J A (i) 



for each A € A. Then 



t—>oo t 



7 := lim - log I V h(t) ] = max 7a . 
t^-oo t \ ^-^ J AeA 



Proof. This follows from a straightforward comparison of the growth of exponentials. □ 
Lemma 3.6. Consider I\ in Iil6\) with f(z) = F^f t (z). Then 

7a = X^V (17) 
i=i 

Here 

7 r = H r (z , r ) 
where (as in the statement of Theorem \l .1^ 

Hr (z) = ± Y ^ + rz- 

and zq^j. is the unique solution to H' r (zQ^ r ) = with ZQ )r G (0, oo). 

Proof. First observe that we can take the contours of integration for Wj in I\ to be a large circle 
containing {0, — 1, . . . , — Aj + 1}. This is because before having applied the identity in Lemma [3.4l we 
could take the Zj contours to be large enough nested circles so that zj, contains 0, — 1, . . . , — Aj + 1. 

Next we observe that the integrals defining I\ match the form of (|22p in Lemma [5.11 with i = £(X) 
and (using io's instead of z's) 

n <(A) 




g(wi,...,w e ) = — —- : det 

mi!m2' • • • 



1 



Aj + Wi - Wj 



i,3 =1 



and Gj (wj ) = H\ j (wj ) . If we can show that the four hypotheses of Lemma 15.11 apply, then the 
result claimed in the present lemma follows immediately. 

By convexity, Gj(wj) has exactly one critical point along Wj G (0, oo). Call this point Wqj. 
Without changing the value of the integrals, we can freely deform the contour of integration for Wj 
to a contour Tj which is defined (see Figure [3] for an illustration) as the union of a long vertical 
line segment going through wqj and a semi-circle enclosing {0, — 1, . . . , — Aj + 1}, with radius large 
enough so that for all r € {0, —1, ... , — Aj + 1} and all Wj G Tj \ {wqj}, \wj — r\ > \wqj — r\ . For this 
choice of contour it is clear that all of the hypotheses of Lemma 15.11 hold. In particular hypothesis 
(2a) holds since Re(wj) is constant along the vertical portion of the contour and decreasing on the 
circular part; and Re[log(w;j + i)] > Re[log(u;o,j + i)] for all Wj / Woj along Tj. □ 

We can now complete the proof of Theorem 11.71 Observe that by combining equation (]15p with 
Lemmas 13.41 and 13.61 we find that 

7fc(l;i/) = max7 A (18) 

Ahfe 
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Figure 3. The contour r,- is a vertical line segment with real part woj joined 
with a semi-circle which encloses {0, —1, . . . , —A,- + 1} as well has the property that 
r G {0, —1, . . . , — Xj + 1} and all Wj G Tj \ {woj}, \wj — r\ > \woj — r\. In words, 
this means that the distance from the points on Tj to the various elements of the set 
{0, —1, . . . , — Xj + 1} is minimized for wqj and strictly larger otherwise. 



where 7a is defined above in (]17p . We claim now that this maximum is attained for all k at 7& and 
hence 7fe(l; v) = jk- If we can show this, then the theorem is proved. 

For k = 1, this result is immediate since there is only one partition of 1. For k = 2 there are two 
partitions to consider A = (1, 1) and A = (2). We claim that for all v > 0, 271 — 72 > 0. This can 
be proved by explicit evaluation. Observe that 

71 = — 1 + v + v\og u, 72 = — 2 + v + \J\ + v 2 — v\og {^y{y + \/l + v 2 )^ . 

The difference f(v) := 72 — 271 goes to zero as v goes to infinity and has derivative log(2i/) — log(^ + 
VT+ v 2 ) which is negative for all v > 0. This shows that f(u) > for all v > and hence 

7 2 (1;^) = max 7 A = 72. 

Al - 2 

Note that we have now shown that 7i(l;z^) < 72(1; ^)/2 and hence we can apply Lemma 1 1.11 to 
show the intermittency of all of the moment Lyapunov exponents. We now proceed by induction on 
k. Assume that for all j < k, we have proved that 7fc(l; v) = 7fe. As a base case we have k = 1 and 
2. By intermittency we know that for any partition of k + 1 aside from A = (k + 1), 7fc+i(l; v) must 
strictly exceed 7A;(1; v) which, by induction, equals J2i 7a ; a s well. By (fT8|) this implies that the 
maximum over A h k + 1 must be attained for the partition A = (k + 1) and hence jk+i (1; v) = 7^+1 
as desired to prove the inductive step and complete the proof of Theorem 11.71 □ 
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3.1. The replica trick. The replica trick is an idea which goes back to Kac [TT] and which has 
received a great deal of attention within the statistical physics community. In its most basic form, 
one hopes to extract the almost sure Lyapunov exponent from the knowledge of all of the moment 
Lyapunov exponents. The reader should be warned that was follows is extremely non-rigorous. 
However, we include it since it is a validation of the replica trick in the context of the one-sided 
parabolic Anderson model (see [10] for this approach implemented in the continuous model discussed 
below in Section H]). 

We would like to compute the almost sure Lyapunov exponent 

71 (1; I/) := lim ^Eflog Z ± (t, vt)}. 

Note that even though we have taken the expectation of log Z\(t, vt), this should not affect the 
value of the almost sure exponent. Recall that for z£C \ R_, 

z k - 1 

logz = lim — - — . (19) 

fe->o k 

We have shown in Theorem 11.71 that 

E[Z 1 (t,ut) k ] e* 7 ^ 1 ^ = e tHk{z °- k) 

for 

fc(fc-3) ( k Y^ t \ k{k-3) T(z + k) 



H k {z) = ^-^ + kz-vlo g [\[{z + i) \=!^L-2l + kz-u] 0t 



a=0 



T(z) 



and zo t k is the unique minimum of Hk(z) for z E (0, oo). This second expression has a clear analytic 
extension in k. 

By using (]19p and interchanging the two limits (without justification) we have 

I e *#fc(zo,fc) _ i 

7i(l:i/) = lim — lim . 

Notice that for k near 0, e tHk( - z °^ w 1 + tHk(zo t k), hence 

~ ( , x ,. 1 ,. tH k {z , k ) 
7i(l;^) = lim — lim 



t-)-oo t fe->o k 

The limit in t can now be taken (since the factors of t cancel) and the limit in k is achieved via 
L'Hopital's rule: 

HmEM. = _3 +z _ 1/ ^) 

where ^{z) = [logr]'(z) is the digamma function. This limit should be evaluated at the unique 
infimum over z G (0, oo) and hence 

li{^) = -- + \^Xz-v^{z)). 



This non-rigorous calculation does yield the proved value (cf. [13] and [12]). 
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4. The space-time continuum parabolic Anderson model 



Consider the solution Zn 
delta function initial data: 



to the multiplicative stochastic heat equation with 



d 



-AZ, 

2 A 



+ PWZ/3, Zp(0,x) = 5 i 



z=Cb 



(20) 



where A is the Laplacian on R, and 5 x= q is the Dirac delta function. The solution Zp(t,x) can be 
thought of as the partition function for a space-time continuous directed polymer in a white-noise 
environment pQ. In fact, under a particular scaling, the parabolic Anderson model considered in 
the previous sections, converges to the SHE [14]. The following formula for joint moments can be 
found by applying this limit transition to Theorem 11.61 For all k > 1 and all x\ < X2 < • • • < x& in 



E 



U 2 ^ 



Xi 



1 



n 



n 



Z- -\~XiZq 



dzi 



(21) 



(2TTi) k J J 11 z a - z b - (3 2 

V ' J J l<a<b<k a ' _x 

where Zj € (3 2 ctj + iM. and a.\ > ct2 + 1 > 0:3 + 2 > • • • > + k — 1. 

The formula of Theorem 11.61 solved the system of ODEs in Proposition 13.11 Likewise, (|2ip solves 
a PDE which is called the 'quantum delta Bose gas' (see Section 6.2 of [3] for a discussion on this, 
as well as remarks on certain gaps in a rigorous statement to this effect). 

From (|2ip it is possible to compute the moment Lyapunov exponents for the SHE (we restrict 
attention now to xi = and to /3 = 1 since general Xi = x and f3 can be achieved from the resulting 
formula via Brownian scaling). The result is that 



Ik 



lim — log 

t— >oo t 



E 



Z(t,0) k 



k 3 -k 
12 



This reproduces Kardar's formula [TU] and agrees with Bertini and Cancrini's Theorem 2.6 and 
remark after it, see [2]. 

This calculation is done by deforming all contours to the imaginary axis and considering the 
growth in t of the various residue terms. As in Theorem 11.71 the Lyapunov exponent comes from 
the (ground-state) term when z± = Z2 + 1 = z$ + 2 = ■ ■ ■ = Zfc + k — 1. There remains only one 
free variable of integration in this residue term, and the main part of the integrand is the following 
exponential: 



cxp 



exp 



+ {z + k-l) 





Deforming the z-contour to tM — ^p- shows that this term behaves like 



exp <{ - ( 1 + 2' + 



lf-k 



{k-lf 



exp < t 



k 3 -k 
12 



from which the result readily follows. 



5. Appendix: Laplace's method 



We following lemma is a version of Laplace's method for computing asymptotics of integrals. The 
proof is an easy modification of the usual proof of Laplace's method [7]. 
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I(t) = jf ••• jf g(z!,...,zt)exp (tY^Gj{zj)\ dzi---dz t . (22) 

Assume that 

(1) For eac/i j, Tj is a closed piecewise smooth contour. 

(2) For each j, there exists z® S Tj such that 

(a) Re[Gj(z)] < Rc[Gj(zQj)] for all z G Tj not equal to Zqj; 

(b) G'j(zQj) = and in a neighborhood of zoj, Gj(z) = Gj(zoj) + c(z — Zoj) r + o((z — zo,j) r ) 
for some r > 2. 

(3) There exists a (non-identically zero) rational function R(z\, . . . , zg) such that in a neighbor- 
hood of(z 0>1 ,...,z j), 

\R(zi, ...,ze)\< \g{zi, • • • ,z e )\. 

(4) There exists a positive constant C such that for all Zj £ Tj (j = 1, . . . , I), 

\g(zi,...,Zi)\ < C. 

Then 

1 1 
Km -logJ(t) = ^G,(^ ). 

°° 5=1 
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